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Large-scale magnetic fields in galaxies are thought to be generated by a turbulent 
dynamo. However the same turbulence also leads to a small-scale dynamo which gen- 
erates magnetic noise at a more rapid rate. The efficiency of the large-scale dynamo 
depends on how this noise saturates. We examine this issue taking into account am- 
bipolar drift, which obtains in a galaxy with significant neutral gas. We argue that, (1) 
the small-scale dynamo generated field does not fill the volume, but is concentrated 
into intermittent rope like structures. The flux ropes are curved on the turbulent eddy 
scales. Their thickness is set by the diffusive scale determined by the effective ambipo- 
lar diffusion; (2) For a largely neutral galactic gas, the small-scale dynamo saturates, 
due to inefficient random stretching, when the peak field in a flux rope has grown to 
a few times the equipartition value; (3) The average energy density in the saturated 
small-scale field is sub equipartition, since it does not fill the volume; (4) Such fields 
neither drain significant energy from the turbulence nor convert eddy motion of the 
turbulence on the outer scale into wavelike motion. The diffusive effects needed for 
the large-scale dynamo operation are then preserved until the large-scale field itself 
grows to near equipartition levels. 

Key words: Magnetic fields; turbulence; Galaxies magnetic fields; ISM:magnetic 
fields; Cosmology: miscellaneous 



1 INTRODUCTION 

Tlie origin of ordered, large-scale galactic magnetic fields 
remains a challenging problem. Magnetic fields in galaxies 
have strengths of order 10~®G, and are coherent on scales 
of several kpc (Beck et al. 1996). These fields can arise, in 
principle, by dynamo generation, from a weak, but nonzero, 
seed field ~ 10~^® — W~^^G, if the galactic dynamo can op- 
erate efficiently enough to exponentiate the field by a factor 
~ 30 - 40 (cf Moffat 1978, Parker 1979, Zeldovich et al. 
1983). However the capacity of presently known turbulent 
dynamo mechanisms to produce the observed galactic fields 
has been debated (Cattaneo and Vainshtein 1991, Vainshtein 
and Rosner 1991, Kulsrud and Anderson 1992, Brandenberg 
1994). 

It has been argued that magnetic noise in the form 
of small-scale fields, builds up much faster than the mean 
field in a turbulent flow. Magnetic noise can result from 
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the tangling of the large-scale field by the turbulence, 
or the operation of a small-scale dynamo. The dominant 
source, when one starts from weak large-scale seed fields, 
is the operation of a small-scale dynamo. Turbulence with a 
large enough magnetic Reynolds number (MRN) , even when 
mirror-symmetric on average, generically leads to an expo- 
nential growth of fields correlated on the turbulent eddy 
scales, independent of the large-scale field. This growth oc- 
curs on the the turbulent eddy turnover timescale which is 
much smaller than the time scale for the growth of the mean 
field. Therefore the kinematic dynamo paradigm will become 
invalid long before the large-scale field has grown anywhere 
near the observed levels. Kulsrud and Anderson (KA) (1992) 
reach the somewhat bold conclusion that the galactic field 
must be of primordial origin! Although there are a number 
of mechanisms to produce small seed magnetic fields (Rees 
1987, 1994, Ratra 1992, Subramanian 1995, Subramanian, 
Narasimha and Chitre, 1995 and references therein ), there 
is as yet no compelling mechanism to produce a field, any- 
where near that required by the primeval hypothesis. 

This problem may disappear if the small-scale dynamo 
generated fields can saturate, due to non-linear back reac- 
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tion effects, in a manner which does not destroy the ability 
of turbulent motions to amplify the large-scale field. We ex- 
amine whether this can indeed happen. Ideally one has to 
consider the MHD dynamo, where both the induction equa- 
tion for the magnetic field, and the Euler equation for the 
velocity field, are solved simultaneously in a self consistent 
fashion. However this is a formidable task at present, even 
for numerical simulation. In this paper we take a first look 
at this difficult non-linear problem, in a simpler fashion by 
isolating the ingredients needed for the small-scale dynamo 
action and investigating under what conditions these get 
suppressed. 

Galaxies have a significant neutral gaseous component. 
As magnetic fields grow, the Lorentz force on the charged 
component will cause a slippage between it and the neutrals. 
Its magnitude is determined by the friction between the com- 
ponents due to ion- neutral collisions. This drift, called am- 
bipolar drift (Mestel and Spitzer 1956, Spitzer 1978, Draine 
1980, 1986, Zweibel 1986), is one important non-linear feed- 
back on both the small and large-scale dynamos, by the 
generated magnetic field. In another paper (Subramanian 
1997 ; Paper II, in preparation), we give a derivation of 
the equations for both mean field and the magnetic correla- 
tions, incorporating the effects of ambipolar drift. We also 
give there the solution of these equations in several con- 
texts. Some pertinent results of this work will be quoted 
here, where needed. We will see below, that the presence of 
neutrals, may be very important, in leading to a saturated 
state for the small-scale dynamo, which preserves large-scale 
dynamo action. 

In section 2, we begin by introducing the turbulent 
galactic dynamo. In section 3, we summarise the properties 
of the kinematic small-scale dynamo action in Kolmogorov 
type turbulence and point out the problem they raise. The 
influence of ambipolar drift on the small-scale dynamo is 
considered in section 4. Section 5 considers the the back re- 
action effects and saturation of the small-scale dynamo, due 
to the Lorentz forces acting on the fluid as a whole. The last 
section contains a discussion and summary of our results. 
We argue that, in galaxies, magnetic noise generated by the 
small-scale dynamo, may indeed saturate in a fashion which 
preserves large scale dynamo action. 



2 THE TURBULENT GALACTIC DYNAMO 

Spiral galaxies are differentially rotating systems. The mag- 
netic flux is to a large extent frozen into the fiuid and so 
any radial component of the magnetic field will be efficiently 
wound up and amplified to produce a toroidal component 
of the field. But this only results in a linear amplification of 
the field: to obtain the observed galactic fields starting from 
small seed fields we must find a way to generate the radial 
components of the field in the galaxy from the toroidal one. 
If this can be done, the field can grow exponentially and one 
has a dynamo. 

The standard picture involves the effects of cyclonic tur- 
bulence in the galactic gas. The galactic interstellar medium 
is assumed to be turbulent, due for example to the ef- 
fects of supernovae randomly exploding in different regions. 
In a rotating, stratified (in density and pressure) medium, 
like a disk galaxy, such turbulence becomes cyclonic and 



acquires a net helicity. Isotropic and homogeneous turbu- 
lence with helicity, in the presence of a large-scale magnetic 
field. Bo, leads to an extra electromotive force of the form 
E = QfBo — 77tV X Bo, where a depends on the helical part 
of the turbulence and rjt is the turbulent diffusion which de- 
pends on the non helical part of the turbulent velocity corre- 
lation function (Krause & Radler 1980, Moffat 1978, Parker 
1979). It should be noted that both the alpha effect, and 
turbulent diffusion, depend crucially on the diffusive (ran- 
dom walk) property of fiuid motion (cf. Field 1996). So, if 
for some reason (see below) the fiuid motion becomes wave- 
like, then the alpha effect and turbulent diffusion will be 
suppressed (Chandran 1996). 

The induction equation, with the extra turbulent com- 
ponent of the electric field, a prescribed large-scale velocity 
field, can have exponentially growing solutions for the large- 
scale field. These have been studied extensively in the lit- 
erature (cf.Ruzmaikin, Shukurov & Sokoloff 1988, Mestel & 
Subramanian 1991, Beck et al. 1996 for a recent review). It 
had been assumed in most earlier works that the turbulent 
velocities do not get affected by Lorentz forces - until the 
mean large-scale field builds up sufficiently. However, this 
does not turn out to be valid due to the more rapid build up 
of magnetic noise compared to the value of the mean field, 
a problem to which we now turn. 



3 THE KINEMATIC SMALL-SCALE DYNAMO 
AND THE PROBLEM OF MAGNETIC 
NOISE 

Split up the magnetic field, B = Bq + (5B, into a mean field 
Bo and a fluctuating component 5B. Here, the mean fleld. 
Bo B >, is defined either as a spatial average over scales 
larger than the turbulent eddy scales or, more correctly, as 
an ensemble average. Kazantsev (1968) was the first to show 
that, even purely mirror-symmetric turbulence, leads to dy- 
namo amplification of the small-scale fiuctuating fields, for 
a sufficiently large magnetic Reynolds number . Some of the 
subsequent work is summarised by Zeldovich et al. (1983). 
The statistical properties of the small-scale field are most 
clearly expressed in terms of the magnetic correlation func- 
tion Mij{r,t) =< (5B'(x,i)5B^(y,t) >, where r = |x - y|. 
Small scale dynamo action is described by the evolution of 
the longitudinal component M_L(r,t) = rV-'Mij/r^, where 
r' = a;' - y\ 

In the case when the turbulent velocity, say vt, has a 
delta function correlation in time (Markovian), it is rela- 
tively straightforward to derive the evolution equation for 
Mij (cf. Kazantsev 1968, Vainshtein & Kichatinov 1986, 
Paper II). Suppose we also assume vt to be an isotropic, 
homogeneous, Gaussian random velocity field with zero 
mean. Specify its two point correlation function by < 
i;^(x,t)tj^(y,s) >= T'^{r)5{t - s), with 

r^(r) = T^^r - (^)] + TllC-^) + Ce^jfT^. (1) 

Here, Tll(t) and TNN{r) are the longitudinal and transverse 
correlation functions for the velocity field and C{r) repre- 
sents the helical part of the velocity correlations, (cf. Landau 
and Lifshitz 1987). If vt is assumed to be divergence- free, 
then, 
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The evolution of Ml is given by 
where we have defined 
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The term involving k in equation (|3|) represents the ef- 
fects of diffusion. The diffusion coefficient includes the ef- 
fects of microscopic diffusion [rf) and a scale-dependent tur- 
bulent diffusion, {Tll{Q) ~ Tll(t)). The term proportional 
to G{r), allows for the rapid generation of magnetic fluctua- 
tions, through shearing action and the existence of a small- 
scale dynamo independent of the large-scale field. [ In the 
evolution equation for Ml given above, we have neglected 
its coupling to the helical part of the magnetic correlation 
due to a scale-dependent alpha effect. In the galactic con- 
text this has a negligible effect. We have also neglected the 
subdominant coupling to the large-scale field (see Paper II) .] 
Suppose V and L are the velocity and correlation 
lengths of the dominant energy carrying eddies of the 
turbulence, which is assumed to have a Kolmogorov en- 
ergy spectrum. For numerical estimates we generally take 

V = lOfcms"^ and L — lOOpc. For Kolmogorov turbu- 
lence, the eddy velocity at any scale I, is w; oc V'^^ , in 
the inertial range. The turbulence is cut off at a scale, say 
la ~ LRe where Re = VL/v is the fluid Reynolds num- 
ber and V is the kinematic viscosity. We will be considering 
a largely neutral galactic gas and for this u is dominated by 
the neutral contribution. We take the neutral-neutral col- 
lision to be dominated by H-H coUisons with a cross sec- 
tion an-H ~ lQ~^^cm~^, leading to a kinematic viscosity 

V ^ Vth{i/nHO'H-H)- For a thermal velocity vth ~ lOfcms"^ 
and a neutral hydrogen number density uh ~ Icm~'^, we 
have u ~ 10^^ cm^s~^, so 

7?, = ^ « 3 X loVioLioo (5) 

where Vio = (V/lOkms"^) and Lioo = (L/lOOpc). The mag- 
netic Reynolds number at the outer scale of the turbulence 
is defined to be Rm ~ (VL/rj). For the Spitzer value of the 
resistivity rj = lO^'CT/lO" iC)"^^^cm^s"\ with turbulence 
parameters as above, Rm = 3 x 10^®. Since vi oc l^^^, the 
magnetic Reynolds number associated with eddies of size / 
scales as Rm{l) = vil/rj = Rm{l/ L)'^^^ ■ 

The evolution of Ml {r, t) has been studied in detail 
by several authors (Zeldovich et al 1983 and references 
therein) for the case of when TLLif) has a single scale. We 
also study in Paper II, using WKBJ technique, the evolu- 
tion of Mi(r, t) for a model Kolmogorov type turbulence. 
We look at the properties of exponentially growing modes, 
ML{r,t) = M„(r)e'""', which are regular at the origin, and 
tend to zero as r ^ oo. Some results of this work and the 
earlier works of Kazantsev (1968) and Zeldovich et al. (1983) 
pertinent to the present context, are summarised here. 

• There is a critical value for the magnetic Reynolds num- 
ber (MRN): Rm = Rc ~ 60, for the excitation of the small- 
scale dynamo. Above this critical MRN the small-scale dy- 



namo can lead to an exponential growth of the fluctuating 
fleld correlated on a scale L. Further, the equations deter- 
mining Rc are the same if we replace {L, Rm) by {I, Rm{l)). 
Therefore, the critical MRN for excitation of a mode con- 
centrated around r ~ Hs also Rmif) — Rc, as expected from 
the scale invariance in the inertial range. 

In the galactic context Rm » Rc', in fact, one also has 
Rmifc) ~ Vclc/ri = Rm/ Re » 1. (Here Vc is the eddy veloc- 
ity at the cut-off scale). Hence, small-scale dynamo action 
excites modes correlated on all scales from the cut-off scale 
Ic to the external scale L of the turbulence. 

• Due to small-scale dynamo action, the fluctuating fleld, 
tangled on a scale I, grows exponentially on the correspond- 
ing eddy turnover time scale, with a growth rate P; ~ Vi/l. 
Since vi ~ V'^^ , the growth rate is P; oc r2/^ and so 
increases with decreasing I. In the galactic context, with 
Rmifc) = Rm/Re >> Rc, the small-scale flelds tangled at 
the cutoff scale grow more rapidly than any of the large-scale 
modes. 

• The WKBJ analysis gives a growth rate Pc = 
{vc/lc)[5/4: - CQ{\n{Rm/Rc))-^] with cq = /12 for the 
fastest mode. Also the growth rate for the small-scale dy- 
namo is only weakly (logarithmically) dependent on Rm, 
provided Rm is large enough. 

• For the parameters adopted above, we have P^^ ~ 
L/V ~ lO'^yr. For the modes tangled at the cut-off scale, 
this time is P~^ ~ 10'*yr. These times are much smaller 
than the time scale for the growth of the large-scale fleld. 

The spatial structure of the dynamo-generated small- 
scale fleld is important in determining how the small-scale 
dynamo saturates. To examine the spatial structure for vari- 
ous eigenmodes of the small-scale dynamo, it is more instruc- 
tive to consider the function w(r,t) =< (5B(x, t).(5B(y, t) >, 
which measures the correlated dot product of the fluctating 
fleld (w(0) =< > ). Firstly there is a general constraint 
that can be placed on wir). Since the fluctuating field is 
divergence free, we have (Kleeorin et al. 1986, Paper II), 



w{r)r dr - 



(6) 
(7) 



since Ml is regular at the origin and vanishes faster than 
r~^ as r — > oo. Therefore the curve r^w{r) should have 
zero area under it. Since w{0) =< (SJi)^ >, w is positive 
near the origin. And the fluctuating fleld points in the same 
direction for small separation. As one goes to larger values 
of r, there must then values of r, say r ^ d, where wir) 
becomes negative. For such values of r, the field at the origin 
and at a separation d are pointing in opposite directions on 
the average. This can be intepreted as indicating that the 
fleld lines, on the average are curved on the scale d. 

• In the case Rm/Re » 1, w{r-) is strongly peaked 
within a region r = rd ~ lc{Rm/Re)^^^^ about the ori- 
gin, for all the modes. Note that r^ is the diffusive scale 
satisfying the condition ri/r^ ~ Vc/lc. 

• For the most rapidly growing mode, w(r) changes sign 
accross r ^ Ic and rapidly decays with increasing r/lc. 
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interpretation as a 
flux rope 




Fig. la 



Fig. lb 



Figure 1. (a) Schematic representation of the the magnetic correlation function ui(r) for the fastest growing mode in the kinematic 
regime. The fact that ui(r) is positive at small r and concentrated into a region of scale near the origin, implies that the field points in 
the same direction on the average for separations smaller than the diffusive scale r^. The anti-correlation tail, with w{r) < 0, is related 
to the vanishing of field divergence; that for points separated by the eddy scale Ic, the field points in opposite directions, on average, (b) 
The pictorial intepretation of ui(r) as a flux rope, with as the rope thickness and Ic as curvature scale of the flux rope. 



• For slower growing modes, with F; ~ vi/l, w{r) extends 
up to r ~ Z after which it decays exponentiaUy. 

We should point out that a detailed analysis of the 
eigenfunctions, for the simple case when the longitudinal 
velocity correlation function has only a single scale, can be 
found in Kleeorin et. al. (1986). Their analysis is also ap- 
plicable to the mode near the cut-off scale in Kolmogorov 
type turbulence. Further, these authors elaborate on a pic- 
torial intepretation of the correlation function, in terms of 
the Zeldovich rope-dynamo (cf. Zeldovich et al. 1983). We 
have shown, schematically, w{r) for the the fastest growing 
mode in Figure la. and its pictorial intepretation interms of 
a flux rope in Figure lb. For the fastest growing mode there 
is only one node where w{r) — 0. For higher order modes 
with smaller growth rates, several nodes for w{r) will occur 
and can be intepreted interms of several scales on which the 
field is curved (cf. Kleeorin et al. 1986, see also Ruzmaikin 
et al. 1989 ). The extent of w{r), (after which it decays ex- 
ponentially), gives the largest scale on which the flux rope 
is curved. 

• If one adopts this intepretation of the w{r) for the var- 
ious modes, the small-scale field can be thought of as being 
concentrated in rope like structures with thickness of order 
the diffusive scale « I and curved on a scale up to ~ Z 
for modes extending up to r ~ Z. 

To end this section we give a qualitative picture of 
the mechanism for the dynamo growth of small-scale fields. 
When the field starts from an arbitrary initial configuration, 
in the kinematic stage, the initial field growth, is just due 
to random stretching by the turbulence, together with flux 



freezing. Eddies of scale I stretch and tangle the field on cor- 
responding scales. The field grows at this stage because the 
stretching of a "fiux rope" leads to a decrease in its cross 
section (due to near incompressibility of the fiow), and hence 
from flux freezing an increase in the field strength. As the 
field grows, the magnetic fleld gets concentrated into smaller 
and smaller scales until diffusion comes into play. 

Consider to begin with the effect of eddies at the cut- 
off scale, Ic, of the turbulence. The initial ampliflcation due 
to purely stretching of the field by these eddies stops when 
the field has been concentrated into a small enough scale Td, 
such that the rate of diffusion across becomes comparable 
to the stretching rate by these eddies. This is when 



(8) 



This gives Vd ~ lc/Rrn'^{lc). We will refer to rd as the thick- 
ness of the fiux rope curved on a scale Ic- Further growth of 
the field can only be achieved by the operation of the small- 
scale dynamo, which exponentiates the field for a MRN 
above Rc. This is explicitly demonstrated by the solutions 
for the kinematic small-scale dynamo discussed above and 
can be thought of as the operation of the Zeldovich - Stretch- 
Twist-Fold - rope dynamo at random locations. 

For galactic gas with Rm/Re >> 1, even the eddies 
at the cut-off scale have an MRN greater than the critical 
value needed for dynamo action. These eddies exponentiate 
the field at a growth rate Vc/lc- Larger eddies of scale I > Ic 
also lead to stretching, twisting and folding of the field at 
a slower rate, vi/l. This leads to dynamo growth of fields 
tangled at scale I, with a slower growth rate Vi/l. The field 
curved on a scale I is also, at the kinematic stage, chopped up 
further by smaller scale eddies (a scale dependent turbulent 
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diffusion ) until its energy can be dissipated by microscopic 
diffusion at scale rd. So, on any flux rope of length I, one 
has smaller scale wiggles until the diffusive scale Vd- 

We emphasise that the time scale for mean-field growth 
is ~ lO^yrs, of order a few rotation time scales of the disk, 
and is much larger than the time scale for the growth of 
the fluctuating fleld ( F^^ ~ lO^yr). Hence, the operation of 
the small-scale dynamo will imply that the magnetic field is 
rapidly dominated by the fiuctuating component. Chandran 
(1996) has shown that the presence of small-scale magnetic 
fields could change hydrodynamic turbulence into magneto- 
elastic waves, with a phase velocity va, where v\ is two 
thirds of the magnetic energy per unit mass. So, if the en- 
ergy density in the small-scale magnetic noise builds up to 
equipartition levels, the fluid motions could become predom- 
inantly wavelike, with a wave period of order the eddy turn- 
over time. This could then lead to a reduced alpha effect and 
turbulent diffusion. Unless the build up of magnetic noise is 
curbed in a way which leaves the turbulence still having a 
diffusive property, large-scale dynamo action will be severely 
affected. We now turn to the effect of ambipolar drift and 
the possible ways in which the small-scale dynamo saturates. 



4 THE EFFECT OF AMBIPOLAR DRIFT 

In a partially ionised medium the magnetic field evolution 
is governed by the induction equation 



(dB/dt) = V X (v, X B - r;V X B), 



(9) 



where Vi the velocity of the ionic component of the fiuid. 
The ions experience the Lorentz force due to the magnetic 
field. This will cause them to drift with respect to the neutral 
component of the fluid. If the ion-neutral collisions are rapid 
enough, one can assume that the Lorentz force on the ions 
is balanced by their friction with the neutrals. Under this 
approximation, the Euler equation for the ions reduces to : 

, _ (V X B) X B 



Pil'ini'Vi 



An 



(10) 



where pi is the mass density of ions, l^in the ion-neutral col- 
lision frequency and v,i the velocity of the neutral particles. 

The ion-neutral elastic scattering frequency is given by 
Vin = Pn < av > j (nii + m„), where p„ is the neutral 
fluid density, and mi,mn are the ion and neutral particle 
masses (cf. Mestel and Spitzer 1956, Draine 1980, 1986). 
We will assume that the galaxy had very nearly primor- 
dial composition in its early stage of evolution: then the 
ions are mostly just protons and the neutrals are mostly 
hydrogen atoms. Elastic scattering occurs with the ion po- 
larising the neutral atom, and interacting with the induced 
dipole. For H — H'^ interactions, in addition to elastic scat- 
tering, there can also be charge exchange reactions, which 
increase the ion-neutral cross-section. Draine (1980) adopts 
the maximum of these two rates, and gives a momentum 
transfer rate coefflcient of < av >~ 3.2 x 10~^cm^s~^ for 
V < 2fcms"\ and < av >x 2.0 x 10~^ (v/kms-^)" ''^ cm^s''^ 
for 2kms~^ < v < lOOOfcms"^. In the galactic disk, we ex- 
pect the gas to have a temperature T < 10* AT with v ~ 
lOfcms"^ and so < av > {H - H+) ~ 10.74 x 10"^cm^s"^ 
Further, the interaction with helium atoms will not give a 
significant addition to the collision rate because the polar- 
isability of helium is less due to its symmetry and helium 



is 4 times heavier than hydrogen. However, helium will con- 
tribute about 25% of the total density of the fiuid. Let 
pH be the proton and hydrogen densities and rii = pu+ jmu 
We then have 

Ph+Ph < av >H-H+ 



pi^i?i 



2m f 



riipn < av >eff 



(11) 



with < av >eff^ 4 X iU^'cm^s" 

For the evolution of the small-scale field, we show in 
Paper II, that ambipolar drift adds to the diffusion coeffi- 
cient, K, a term proportional to the energy density in the 
fiuctuating fields. This changes i] to an effective value 



■w(0,t) <5B.SB> 
■ rj + = ri + 



(12) 



(ynpiVin Gnpiiyin 
and replaces k in equation (^ for ML{r,t), by a new 
tjv = Vambi + Tll{Q) — Thhir) ■ Oue can define an effective 
magnetic Reynolds number, incorporating the effect of am- 
bipolar drift, for fiuid motion on any scale of the turbulence 
by 



Ramhi (^) 



vd 

'Hambi 



VllGnpiVin 

< 5B.SB >' 



(13) 



where vi = {l/Ly^^V as before. 

As the energy density in the fiuctuating field increases, 
Rambi(l) decrcases. Firstly, this makes it easier for the field 
energy to reach the diffusive scales, from a general initial 
configuration. After this stage, the initial ampliflcation due 
to purely stretching stops, and further growth of the fleld 
can only be achieved by the operation of the small-scale dy- 
namo. If, as the fleld grows, Rambi{L) decreases sufficiently, 
a stationary state with dML/dt = could, in principle, be 
achieved. In such a state. Ml is independent of time. So, the 
condition on the critical MRN for the stationary state to be 
reached, will be identical to that obtained in the kinematic 
stage. That is, if Rambi{L) decreases to a value Rc ~ 100, 
dynamo action will stop completely. 

However, for galactic turbulence, 

3pm„l _ Q{1) 

/(O ' 



(14) 



f{l) 2pnVl 

where f{l) — Bf / {A-Kpnvf ) is the ratio of the local magnetic 
energy density of a flux rope curved on scale I, to the turbu- 
lent energy density p„W; /2 associated with eddies of scale I. 
Using the value of determined in Eq. ( pT| ) and putting 

in numerical values we get 



Q{1) = ~ 1.8 X 10*n_2(l)^/^LiooVlo' 

2p„vi L 



(15) 



where n_2 = {ni/lQ~^cm~'^), and we have assumed a Kol- 
mogorov scaling for the turbulent velocity fluctuations. 

One can see from Eq. (^) - (^) that, for typical pa- 
rameters associated with galactic turbulence, the MRN in- 
corporating ambipolar drift is likely to remain much larger 
than Rc for most scales of the turbulence, even when the 
fleld energy density becomes comparable to the equiparti- 
tion value. So ambipolar drift by itself cannot saturate the 
small scale dynamo. Rather, one expects the fleld to con- 
tinue to grow rapidly, even taking into account ambipolar 
drift. Note also that the growth rates for the small-scale dy- 
namo generally depends only weakly on the MRN, provided 
the MRN is much larger than Rc- (see section 3, Kleeorin et 
al. 1986, Ruzmaikin et al. 1989). Therefore, we still expect 
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the small-scale dynamo-generated field to grow almost ex- 
ponentially on the eddy turn around time scale, as long as 
R » Rc. 

The spatial structure of the fluctuating field will also re- 
main ropy, as argued in section 3, as long as RamU » Rc- 
However, as the field strength in a fiux rope, curved on a 
scale > grows to near equipartition with the turbulent 
energy associated with eddies smaller than I, these smaller 
scale eddies would no longer be effective in causing "turbu- 
lent diffusion" of the larger scale field. So the thickness rd{l) 
of a flux rope curved on a scale I is determined eventually 
by demanding that the ambipolar diffusion timescale across 
rd{l) becomes comparable to the stretching timescale l/vi; 
that is, 



rdiO/Varnbi ^ I / Vl . 



(16) 



This determines the thickness to be rd{l) ~ V^ambiW- Since 
Rambi(l) >> 1, wc expect flux ropcs to remain relatively 
thin with a thickness ~ '/^ambiW << ^) even taking ac- 
count of the ambipolar drift. 

In summary, we have argued here that, in galaxies, the 
small-scale dynamo continues to exponentiate the field fluc- 
tuations even in the presence of ambipolar drift. This fluc- 
tuating field however does not fill the volume but is con- 
centrated into intermittent rope like structures. The ropes 
are curved on the turbulent eddy scale and their thickness 
is set by the diffusive scale rd{l) determined by the effective 
ambipolar diffusion. 

We have to consider how other non-linear feedback pro- 
cesses could limit the growth of this small-scale dynamo gen- 
erated field. 



5 SATURATION OF THE SMALL-SCALE 
DYNAMO 

5.1 Inefficient random stretching and damping 

The first of these restraining processes is the reduction in 
the efficiency of stretching of a flux rope as the field in the 
rope, say Bp, grows in strength. A turbulent eddy of scale / 
produces a correlated winding up of the field for a time of 
order l/vi. Note that this is just a consequence of the induc- 
tion V X (v X B) term and velocity shear. No dynamics is 
involved. However, suppose the growing tension component 
of the Lorentz force on the flux rope tangled on this scale, 
can untangle (or straighten) the rope and damp away its 
wrinkle, on a comparable timescale. The random stretching 
of the flux rope by these eddies will be suppressed and the 
small scale dynamo will not operate efficiently. 

To estimate the straightening time-scale, say ts, we 
have to look at the dynamics of the flux rope. The mo- 
tion of the fiux rope through the surrounding medium is 
influenced not only by the tension force but also by fric- 
tion. (We are neglecting, for now, the effects of rotation 
and gravity). Suppose Vs is the velocity associated with the 
untangling or straightening of the field in the flux rope. 
If one were to ignore the effect of friction, we would have 
Vs ^ va ~ {Bp/A-KpnY^^ . However, the effect of drag leads 
to a smaller "terminal" velocity, which can be estimated by 
equating the tension force to the drag on the flux rope. Also, 
the work done by the rope against the frictional drag leads to 



a damping of the energy associated with the wrinkle of the 
flux rope on the straightening timescale, l/vs, which is com- 
parable to the rapid eddy turnover time when Vs ~ We 
illustrate schematically, this process of flux rope straighten- 
ing and damping in Figure 2a. 

The drag force per unit length on the flux rope is 
^ Cdpuvlrd/2 (cf. Parker 1979 Eq. 8.59). The coefficient 
Cd depends on the fluid Reynolds number on the scale of 
the radius of the flux rope, that is on JJ^'"* — VsTd/v- Note 
that the drag formalism is not well developed for the case 
when the external medium is itself turbulent (cf. Parker 1979 
Section 8.7). One expects eddies with a scale smaller than 
the rope radius to have a different effect compared to larger 
eddies. In our case, rd « I, and so one expects momentum 
transfer by eddies with scale < rd to be subdominant com- 
pared to the turbulent drag induced by larger eddies. Here 
we adopt the drag formula given above, with the assumption 
that the complications mentioned above do not drastically 
change the results. 

Equating the magnetic tension component of the 
Lorentz force in the rope and drag then gives 



/ 2\ Cd 2 



(17) 



Here, we have taken the curvature radius associated with 
the field tangled on scale I to be 1/2. The dynamo action 
on these scales will be affected when ts ~ l/vs ~ l/vi or 
when Vs ^ Vt- Since, in this case, the flux rope is able to 
straighten and damp (due to friction) its wrinkle (curvature) 
on a timescale comparable to the stretching timescale l/vi. 
So, the random stretching of the fleld by eddies of scale /, 
will become inefficient when the peak fleld hasgrown to a 
value found by substituting Vs ~ vi in Eq. (|l7|). Straight 
forward rearrangment of the various quantities in Eq. ( 
) then implies an upper limit on the magnetic field in the 
flux rope, of 



d2 

Ml)- " 



CdJ_ 

ipnvf/2) ~ 4-Krd' 



(18) 



Since we expect rd < I, the peak fleld in the rope can exceed 
equipartition values. 

When ambipolar drift is causing the flux rope to 
thicken, we have l/rd = Rai^^iil), which itself depends on 
the magnetic field strength in the rope. Further, the drag 
co-efficient, Cd, depends on i?^"'', so on rd, and therefore 
implicitly on the field in the rope. 



The value of Cd decreases from about 30 at Rl 



1 to 



4 at Rl""^ ~ 10, 1 at Rl"'^ - 10^ and ranges between 1-0.1 
for larger values (cf. Parker 1979 Sec. 8). For our problem, 
an adequate estimate of the drag co-efficient is obtained by 
taking Cd ~ 10/^_R?'"'(rd). Taking Vs ~ w;, we have i?™"* = 

ReH){rd/l)- Using l/rd = Rl^r^biH), and re-arranging, we 
obtain 



3 1 n^'/^r-r^/^'L^/V"^/^ 



(19) 



One can now go back and estimate the Reynolds number 
Rl""^. We get 



Re. 



392{l/Lr/''nZrV,TLlil 



(20) 
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Figure 2. (a) A schematic illustration of a curved flux rope, with a radius of curvature I and thickness r^j. Magnetic tension acts to 
straighten the rope and aerodynamic drag damps the magnetic energy associated with the wrinkle in the rope. This leads to inefHcient 
random stretching when Vs becomes comparable to t); . (b) The collapse of flux loops is illustrated schematically. Such collapse results in 
an irreversible removal of small-scale magnetic noise and limits the value of N. 



At the cut of! scale l^, i?™'' ~ 0.4, while at the outer scale 
L, iZg"'' ~ 392. So the approximation used for the drag co- 
efficient should be reasonably accurate. 

We see from Eq. (^9|) , that when the field in the flux 
rope of scale I grows to a few times the equipartition field 
associated with eddies on that scale, the random stretching 
of the rope will become inefficient, limiting its further growth 
by dynamo action. Since larger eddies carry larger energy, 
this also implies that the the saturated value of the fleld in 
a flux rope, tangled on a scale I will be larger for larger I. 

The average energy density in the fleld tangled on the 
scale I, however need not exceed the energy density of the 
turbulent eddies because of the ropiness of the fleld. Since 
the ropes have a thickness of order and are curved an 
a scale I, the average energy density contributed by a flux 
rope to a sphere of radius I is 



Br, 



Stt (47r/3)P Stt 4 P ' 



(21) 



Note that, although oc (l/ra), and can exceed equipar- 
tition value, the average energy density of a flux rope in a 
correlation volume Eb{1) oc Bp{rd/Vf (x (rd/l)- This can be 
much smaller than equipartition because Vd/l « 1. 

In principle, the length of the flux rope tangled on a 
scale I may be larger than the value ttI assumed above, and 
equal to say Nnl. In this case, the ratio of the average mag- 
netic energy density of the fleld tangled on scale I to the 
turbulent energy on the same scale is 



F{1) = 



Eb{1) 



iV(^)(^). 



Putting in numerical values, we have 



(22) 



F(0 =0.9 X lQ-'^NnZy\Y) 

1j 



/6r-3/4 



1/4 



(23) 



We can also ask if the energy dissipated in ambipo- 
lar drift and frictional drag is comparable to the turbu- 
lent power in the case when the small scale fleld satu- 
rates due to the stretching constraint. Since the ambipo- 
lar drift rate across the rope is comparable to the stretch- 
ing rate, vi/l, the power dissipated by ambipolar drift is 
{riambi/rl{l))EB{l) ~ Eb{1){vi/1). The power lost due to 
frictional drag is also ~ Eb{1)(1/vi). So the ratio of the to- 
tal power dissipated to the turbulent power is 



Pd 
Pt 



= 2 



Eb{1){vi/1) 
{pnvf/2){v,/l) 



2F{1) 



(24) 



Both the energy density in the flux ropes, and the en- 
ergy dissipated due to ambipolar drift and friction depend 
on how large can get, or how many flux ropes of scale I 
are packed into a correlation volume of radius /, in the flnal 
saturated state. 



5.2 The limiting effect due to flux loop collapse 

An important process, which limits A'' from becoming too 
large, enters as the peak field in the flux rope increases to the 
value given by Eq. (p^. As one tries to pack more and more 
flux ropes curved on a scale I into a volume of same scale, 
the probability that the rope intersects itself, or another 
co-habiting rope, increases. This will result in loops of flux 
of radius I. Also, note that the the very process by which 
the small-scale dynamo may operate, viz. via the stretching- 
twisting-folding actions associated with the Zeldovich rope 
dynamo, will generically result in flux loops of scale I. 
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These loops of flux, when not being stretched, can col- 
lapse to a small radius (cf. Deluca, Fisher & Patten 1993, 
Vishniac 1995) .This is schematically illustrated in Figure 
2b. The time-scale for the collapse is of order l/vs, which 
will be ~ l/vi, the eddy turnover time, since Vs ~ vi when 
the peak field in the rope is given by Eq. (p^. The collapse 
of the loop will be halted when it has shrunk to a radius 
comparable to its thickness. At this stage, the effects of dif- 
fusion could convert most of the remaining energy in the 
loop into kinetic energy and heat. This process results in 
the irreversible removal of energy from the small-scale mag- 
netic field tangled on scale I at a rate Vi/l, comparable to 
its rate of build up by turbulent stretching, after the field 
has grown sufficiently. A dynamical equilibrium for N can 
then result, whereby "old" loops of flux, which have peak 
flux given by Eq. (|l8|), collapse and are destroyed at a rate 
vi/l, to be replaced by newly created loops, at the same rate, 
which are just reaching the saturated value of the peak flux. 
This picture of dynamical equilibrium leads one to conjec- 
ture that, in the saturated state, the average length of flux 
ropes tangled on scale /, in a volume of radius I, cannot grow 
much larger than ttZ. That is the value of A'', can not grow 
much larger than unity. 

Given that N is not too large, one can see from Eq. 
( p^ that the power dissipated in ambipolar drift is much 
smaller than the turbulent power. This occurs for all scales 
except near the cut-off scale. In fact, one only needs N 
less than about 100 for the eddies at the energy-carrying 
scale to be left unaffected (undamped) by the growth of the 
small-scale field. Also, from Eq. (^^, one can see that the 
average energy density of the generated small-scale field is 
much smaller than the average energy density in the turbu- 
lence. So, any wave-like motion induced by the presence of 
the small-scale field (cf. Chandran 1996) will have a period 
larger than the eddy turn around time. This implies that 
such tangled small-scale fields do not change the diffusive 
nature of the turbulence. For these reasons, the large-scale 
dynamo can still operate to generate the mean field. 



5.3 The limit on the field in flux ropes due to the 
external pressure of the gas 

There is one caveat to the above discussion. We have as- 
sumed that the field in the rope can grow sufficiently so 
that the tension in the rope begins to play an important 
role in the rope dynamics. However, there is an upper limit 
to the growth of the magnetic field in the ropes, from the 
effect of its magnetic pressure on the dynamo process. Due 
to the increasing importance of this pressure, stretching of 
field lines can lead to a partial decrease in fiuid density in 
the ropes rather than a decrease in the rope-cross section 
and the associated increase in the rope magnetic field (cf. 
Vishniac 1995). An upper limit to the magnetic pressure in 
the ropes is given by the external pressure Pext- This implies 
that the field in the rope, is limited to 

\l/2 



B„ < (87rPe. 



(25) 



Whether the field in the ropes will be limited by inefficient 
random stretching as given by Eq. (^9|) or by external pres- 
sure (Eq. (pi])) will depend on the parameters of the prob- 
lem. The rope field will be given by the lower of the limiting 
fields implied by Eqs. (|l|) and (H). 



The total pressure of the interstellar medium in a galaxy 
could contain a number of components; a thermal compo- 
nent, pressure due to turbulence itself and possibly due to 
non-thermal "cosmic- rays" . The ratio of the gas pressure to 
the turbulent energy density is 



PJEt ~ 1.7(r/10*A-)V^io'- 



(26) 



If Pcxt in a galaxy is a factor F times the gas pressure, then 
the peak field given by Eq. ( |l9| ) begins to exceed that given 
by Eq. (psj), when the ion density exceeds a critical value 
n1. This critical value is given by 



6.4 X 10-^(^^)2/\ 



T 



■104 if 



7/3 



(27) 



For a larger ion density than nl, the peak field for the 
fiux ropes tangled on the largest scale L saturates to a value 
Bp ~ {%-KPea:tf^'^, lower than given by Eq. (|l|) . So, for 
the small-scale dynamo-generated field to saturate to sub- 
equipartition level by the processes described above, the 
galactic gas has to be predominantly neutral with an ion 
density less than nl. Note that the value of nl is critically 
dependent on the turbulence parameters that obtain in the 
ISM of the galaxy, especially the turbulent velocity scale. 
For example for V = 5kms~^ , and all other parameters as 
above, the critical density becomes nl ^ 0.5cm~'^; so even 
for an ionised hydrogen density as large as 10% — 50% of a 
neutral density, taken here to be lcm~^), saturation could 
occur due to inefficient stretching. 

Let us now ask what happens if the ion densities exceed 
the critical value n^. In this case {Bp/Sir) < Pext and the 
stretching constraint can only be satisfied only if the fiux 
rope thickens further than the value implied by ambipolar 
drift, to a radius rd = R ^ l{Cd/4:Tv){pnvf / {2Pext)) ■ If the 
fiux ropes can thicken to this radius, the average energy 
magnetic energy density will be 



Eb{L) 
Et 



N- 



1287rF 



-Mt 



7.3 X IQ-^NCl 



(28) 



where Mt is the Mach number of the turbulence and we 
have adopted the parameters 7^ ~ 4 and Pg/Er ^ 1.7, given 
above for the numerical estimate. The value of N should 
again be limited by the collapse of loops as discussed above. 

However, it is not clear if fiux ropes can thicken further 
than the radius implied by ambipolar drift, when dynamo 
action begins from weak seed fields. One possibility is that 
the magnetic pressure in the ropes acting on the fiuid as 
a whole can thicken the rope. But this can only happen if 
the pressure in the rope becomes larger than the pressure 
outside, at least temporarily. Even if this were possible, as 
the flux rope thickens, fiux freezing leads to a decrease of 
the field strength in the rope, and a consequent decrease in 
the thickening rate. The fluid pressure in the rope will also 
decrease as the rope expands. Note that this problem does 
not arise when ambipolar drift is causing the thickening. In 
a predominantly neutral medium, the ion pressure is much 
smaller than that due to neutrals. And the Lorentz force 
term in the Euler equation for the ions can be much larger 
than the ion pressure gradient term, and cause a relative 
drift of the ions with respect to the neutrals, carrying the 
field, and hence thickening the fiux rope. 
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If the ion density is larger than and the flux ropes 
cannot thicken sufficiently to resist stretching, then it is not 
clear how exactly the small scale field saturates. The field 
may get packed locally into a radius ~ i? by the folding mo- 
tions associated with the turbulence, and then resist further 
stretching. On the other hand, the small-scale field may only 
saturate if the length of flux ropes, N, increases sufficiently 
to achieve equipartition with the turbulence. In either case 
the small-scale field will be highly intermittent. The effects 
of reconnection of this highly intermittent field (cf. Vishniac 
1995, Lazarian and Vishniac 1996) could be important in 
deciding if the turbulence can still lead to large-scale dy- 
namo action. In an interesting paper which came to our 
notice during the completion of the present work. Black- 
man (1996) discusses the possible effects of reconnection in 
greater detail, albiet in the case where the fiux ropes are 
assumed to have a thickness ~ R, and assuming the gas is 
largely ionised. 



6 DISCUSSION AND CONCLUSIONS 

The large-scale galactic field is thought to be generated by 
a turbulent dynamo. However the same turbulence will pro- 
duce magnetic noise at a more rapid rate. We have exam- 
ined whether the Lorentz forces associated with the growing 
small-scale fields, can lead to their saturation, in a manner 
which preserves large-scale dynamo action. In doing this, we 
have also taken account of the ambipolar drift induced by 
the presence of a neutral component of the galactic gas. 

The saturated state of the small-scale dynamo gener- 
ated field, which we have motivated, does indeed preserve 
large-scale dynamo action. The crucial property of the small- 
scale dynamo generated field which allows this to happen is 
its spatial intermittency. The field can build up locally to a 
level which will lead to small-scale dynamo saturation, while 
at the same time having a sub-equipartition average energy 
density. 

Numerical simulations of dynamo action due to mirror- 
symmetric turbulence (Meneguzzi et al. 1981) or convection 
(Brandenburg et al. 1996) have indeed hinted at a saturated 
state of the small-scale dynamo as described above; a mag- 
netic field concentrated into flux ropes, occupying a small 
fraction of the fluid volume, having peak fields comparable 
or in excess of equipartition value but average magnetic en- 
ergy density only about 10% of the kinetic energy density. 

We have described in section 4 and 5 the approach to 
this saturated state. As we noted in section 4, for conditions 
appropriate to galactic gas, the effective magnetic Reynolds 
number, even including ambipolar diffusion, is much larger 
than a critical value needed for small-scale dynamo action. 
However, in such a case, as the the small-scale field grows 
in strength, it continues to be concentrated into thin ropy 
structures, as in the kinematic regime. These fiux ropes are 
curved on the turbulent eddy scales, while their thickness 
is set by the diffusive scale determined by the effective am- 
bipolar diffusion. The growing magnetic tension associated 
with the curved flux ropes, acts to straighten them out. Fric- 
tional drag damps the magnetic energy associated with the 
wrinkle in the rope. Also, small-scale flux loops can collapse 
and disappear. These non-local effects operate on the eddy 
turnover time scale, when the peak field in a flux rope has 



grown to a few times the equipartion value. Their net effect 
is to make the random stretching needed for the small-scale 
dynamo inefficient and hence saturate the small-scale dy- 
namo. However, the average energy density in the saturated 
small-scale field is sub equipartition, since it does not fill 
the volume. Such fields neither drain significant energy from 
the turbulence, nor convert eddy motion of the turbulence 
on the outer scale to wave-like motion. The diffusive effects 
needed for the large-scale dynamo operation are then pre- 
served. This picture of small-scale dynamo saturation ob- 
tains when the ion density is less than a critical value of 
n1 ~ 0.06 - 0.5cm-^(n„/cm-^)^/^ 

For very large ion densities Ui > n1, the small-scale field 
is expected to saturate only when its energy density grows 
comparable to that of the turbulence. This is because, in 
this case the peak field, cannot grow sufficiently, (without its 
pressure exceeding the interstellar pressure) , for the stretch- 
ing constraint to apply. However, the spatial structure of 
the small-scale field, is still likely to be highly intermittent. 
The large-scale dynamo action will depend on how such a 
field responds to turbulent motions, especially whether the 
field can reconnect efficiently (cf. Vishniac 1995, Blackman 
1996). 

We discuss briefiy, the implications of the above re- 
sults for the origin of galactic fields. The viability of the 
saturation mechanism discussed here for limiting magnetic 
noise depends on the ionisation of the gas and the turbu- 
lence parameters. In the context of the Galaxy, a study 
of warm clouds by Spitzer and Fitzpatrick (1993) gives a 
range of electron densities for the clouds, with an average of 
~ 0.07cm~"^. They also deduce an average neutral density 
~ 0.2cm~'^. For such parameters, the magnetic noise will 
indeed saturate if the turbulent velocity is ~ 5kms~^ and 
other parameters are as in Eq. ( p7| ) . 

In the case of a young galaxy with mass ~ 10^^ Mq that 
has just collapsed into a lOkpc sized region, the average den- 
sity is larger, and is ~ lcm~^ . The average column density 
for HI is ~ lO^^cm"^. The damped Ly-a systems seen in the 
spectra of high redshift quasars, and thought to be young 
galaxies do indeed have such HI column densities (cf. Wolfe 
1995). The ionisation fraction is more uncertain. For such a 
high HI column density as inferred above, the gas is expected 
to become self shielded to external ionising flux. However, 
ionisation will result from UV emission from young stars 
embedded in the gas, whose importance depends on the un- 
certain star formation rates and stellar mass functions. From 
an observational point of view, a recent study of metal lines 
in these systems (Lu et al., 1996) deduces an upper limit to 
their electron density, consistent with the average electron 
density in warm clouds in the Milky Way ^ O.OTcm""^ men- 
tioned above. Further velocity widths deduced from 21cm 
absorption due to neutral hydrogen in several damped Ly-a 
systems, also limit a turbulent component to the line width 
to be about lOfcms"'^ (cf. De Bruyn et al. , 1996). So, here 
also, the densities and turbulence parameters are expected 
to be in the range wherein the small-scale dynamo generated 
fields can saturate due to the tension forces, in a way which 
preserves large-scale dynamo action. 

In our analysis, so far, we have ignored the small-scale 
field generated by the tangling of the large-scale field by the 
turbulence. Even when the dynamo generated small scale 
field has saturated, this will provide an additional source 
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of small-scale magnetic noise. As the large-scale field grows, 
so does this component of the small-scale field with an en- 
ergy density ultimately decided by the nature of the MHD 
turbulence (cf. Zeldovich et al. f983). We hope to return 
to this issue in a later work. One then expects two com- 
ponents to the small-scale magnetic field in the interstellar 
medium of a galaxy. First, a ropy, intermittent component, 
with flux ropes curved on scale L ~ f OOpc, say, and thickness 
Td ~ fO~^L ~ fpc, with peak field a few times equipartition. 
Second, a more diffuse small-scale field related in strength 
to the large-scale field, ft would be interesting to search 
for both of these components in the interstellar medium of 
galaxies. 

The effect of ambipolar drift, together with the small- 
scale dynamo, can also influence galactic magnetic field gen- 
eration, indirectly, in another fashion. Note that any dynamo 
needs a seed field to act upon. Since the small-scale dynamo 
acts to generate fields more rapidly than the large-scale dy- 
namo, the magnetic noise so generated may itself provide a 
significant seed for the large-scale dynamo (cf. Beck et al. 
f994). ffowever the small-scale dynamo also leads to highly 
ropy fields with a rope thickness, r^. In a fully ionised gas 

= ~ L/rH'^ will be very small, since Rm » 1- The 
overlap of such a field with a large-scale dynamo eigenfunc- 
tion will be small, ffowever, if small-scale dynamo action 
proceeds in the presence of neutrals, r^j = ~ ^/^tmbi 

in general, since Rambi « Rm- So, when neutrals are 
present, the small-scale dynamo-generated magnetic noise 
will provide a more coherent seed field, for large-scale dy- 
namo action. This will act to shorten the timescale for the 
generation of large-scale galactic magnetic field, to the mi- 
crogauss level, at higher redshift. 

Note that, after recombination, the residual ionisation 
fraction of the intergalactic medium drops to about f — 
fO"^ One may be tempted to apply some of the results ob- 
tained here to the first generation of objects which collapse 
at high redshifts. The limitations of our semi-quantitative 
arguments, and our assumption of a homogeneous galac- 
tic interstellar medium, in reaching the above conclusions, 
needs little emphasising, ft would also be fruitful to find 
ways of incorporating, more fully, the dynamics of the ve- 
locity correlations, as we have done for the magnetic correla- 
tions. This full MHD turbulence problem appears formidable 
at present. Nevertheless, the results obtained here encour- 
age the belief that that the turbulent galactic dynamo could 
indeed be made to produce large-scale fields, in the presence 
of a significant neutral component. 
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